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Aims. Neutron star crust, formed via accretion of matter from a companion in a low-mass X-ray binary (LMXB), has 
an equation of state (EOS) stiffer than that of catalyzed matter. At a given neutron star mass, M, the radius of a star 
with an accreted crust is therefore larger, by AR(M), than for usually considered star built of catalyzed matter. 
Methods. Using a compressible liquid drop model of nuclei, we calculate, within the one-component plasma approxi- 
mation, the EOSs corresponding to different nuclear compositions of ashes of X-ray bursts in LMXB. These EOSs are 
then applied for studying the effect of different formation scenarios on the neutron-star mass-radius relation 
Results. Assuming the SLy EOS for neutron star's liquid core, derived by Douchin & Haensel (2001), we find that at 
M — 1.4 Mq the star with accreted crust has a radius more than 100 m larger that for the crust of catalyzed matter. 
Using smallness of the crust mass compared to M, we derive a formula that relates AR(M) to the difference in the 
crust EOS. This very precise formula gives also analytic dependence of AR on M and R of the reference star built of 
catalyzed matter. The formula is valid for any EOS of the liquid core. Rotation of neutron star makes AR(M) larger. 
We derive an approximate but very precise formula that gives difference in equatorial radii, A.R cq (M), as a function of 
stellar rotation frequency 
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1. Introduction 

Establishing via observations the mass versus radius rela- 
tion for neutron stars is crucial for determination of the 
equation of state (EOS) of dense matter. It is expected that 
finding even a few points of the M(R) curve could severely 
limit the range of considered theoretical models of matter in 
a liquid neutron star core. General method of determining 
the EOS at supranuclear densities (i.e. densities larger than 
the normal nuclear density po = 2.7xl0 14 g cm~ 3 ) from the 
M(R) curve was developed bv lLindblom I ( 1992f ). A simpli- 
fied method of determining a supranuclear segment of the 
EOS from three measured points of the M(R) has been re- 
cently described in (jOzel et al. I . l2010h . In these studies, the 
EOS of neutron star crust is considered as well established 
and fixed. This assumption is valid, within relatively small 
uncertainties, provided the crust matter is in full thermo- 
dynamical equilibrium (catalyzed matter, corresponding at 
T = to the ground state of the matter). Such a condition 
is not fulfilled for the crust formed by accretion of matter 
onto a neutron star surface, f rom a companion in a low-m ass 
X-ray binary (accreted crust . [Haensel fe Zdunik! (EHoiE)). 
The EOS of accreted crust is st i ffer tha n that of the ground- 
state one (|Haensel fc Zdunikl , Il990b| ). Therefore, at the 
same stellar mass. M, the radius of the neutron star with 
an accreted crust is larger by AR(M) than that of the star 
with crust built of catalyzed matter (catalyzed crust). We 
expect t hat the millisecond pulsars, spun-up by ac cretion in 
LMXBs (jBhattacharva fc van den Heuvel 1 . 1 1 9 9 lh have ac- 
creted crusts, different from the catalyzed crusts of pulsars 
born in supernova explosions. 
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In the present note we study the effect of the forma- 
tion scenario on the radius-mass relation for neutron stars. 
Formation scenarios and corresponding equations of state 
of the crust are presented in Sect. [2] In Sect. [3] we cal- 
culate AR(M) for non-rotating stars, assuming different 
crust formation scenarios. Numerical results for AR(M), at 
fixed EOS of the liquid core, are presented in Sect. 13. ll In 
Sect. 13.21 we derive an approximate, but very precise, gen- 
eral formula that relates the difference in radii to the differ- 
ence in the EOSs of the crust. The dependence of AR(M) 
on the EOS of the core enters via a dependence on the 
mass and radius of the reference star with catalyzed crust. 
Rotation of neutron star increases AR(M), and in Sect. [4] 
we study AR(AI) (defined as the difference between the 
equatorial radii) for rotating neutron stars. Section [5] con- 
tains a summary of our results, discussion, and conclusion. 

Calculations that led to the present note were inspired 
by questions from the audience during a talk by one of the 
authors (P.H.) at the CompStar Workshop "Neutron star 
physics and nuclear physics" , held at GANIL, Caen, France, 
February 14-16, 2010. 

2. Scenarios and equations of state 

2.1. Catalyzed crust 

We assume that neutron star was born in a core-collapse 
supernova. Initially, the outer layers of the star are hot 
and fluid. Their composition corresponds to nuclear equi- 
librium, because at T > 10 10 K all reactions are sufficiently 
rapid. The crust solidifies in the process of cooling of a 
newly born neutron star. One assumes, that while cooling 
and solidifying, the outer layers keep the nuclear equilib- 



2 



J.L. Zdunik, P. Haensel: Formation and R(M) for neutron stars 



rium, and after reaching the strong degeneracy limit the 
state of the EOS of the crust is well approximated by that 
of the cold catalyzed matter, corresponding (in the T = 
limit) to the ground of the mat ter (for a detailed discussion, 
see, e.g.. lHaensel et al. I (|2007tl ). This EOS of the crust will 
be denoted as EOS[cat]. 

2.2. Accreted crusts 

We assume that neutron star has been remaining ~ 10 8 — 
10 9 yr in a LMXB. Its crust was formed via accretion of 
matter onto neutron star surface from a companion in the 
binary. Typical accretion rates in LMXBs arc ~ 10~ 10 — 
10~ 9 M Q yr -1 . Therefore, as a result of the accretion stage, 
the original catalyzed crust formed in scenario described 
in Sect. I2.ll (of typical mass ~ 10~ 2 M Q ) has been fully 
replaced by the accreted one. 

Accretion and X-ray bursts. Many LMXBs are sites of 
type I X-ray bursts (hereafter: X-ray bursts) which are ther- 
monuclear ^Jhishej3_jii_thB_j^ 

tron star (IWoosley fc Taam I . Il97rl iMaraschi fc Cavaliere I . 
119771 I Joss I 119771) . In the simplest model of the X-ray 
bursts, accreted matter, composed mostly of 1 H, cumulates 
on the star surface and undergoes compression due to the 
weight of the continuously accreting matter. The accreted 
layer is also heated by the plasma hitting the star surface 
and transforming its kinetic energy into heat. Compressed 
and heated hydrogen layer burns steadily, in its bottom 
shell with p ~ 10 5 g cm~ 3 , into 4 He. The helium produced 
in the burning of hydrogen is accumulating in a growing He 
layer beneath the H-burning shell. After some recurrence 
time (typically ~ hours), the helium burning is ignited at 
the bottom of the He layer, typically at p ~ 10 6 g cm -3 . 
The He burning starts in a strongly degenerate plasma 
(temperature ~ 10 8 K and p ~ 10 6 g cm -3 ). Therefore, He 
burning is thermally unstable and proceeds initially in an 
explosive detonation mode, with local temperature exceed- 
ing 10 9 K, and burns the overlaying He and H layers into 
elements with mass number A ~ 50 — 100 (see next para- 
graph). Finally, the thermonuclear explosion develops into 
a thermonuclear flash of the surface layer, observable as 
an X-ray burst. After ~ minutes H-He envelope has been 
transformed into a layer of nuclear ashes. The energy re- 
leased in thermonuclear burning has been radiated in an 
X-ray burst. Continuing accretion leads again to the cumu- 
lation of the H-He fuel for the next X-ray burst and the 
cycle accretion-burst repeats in a quasi-periodic way every 
few hours. 

Ashes of X-ray burst. The composition of ashes from 
thermonuclear burning of an accreted H-He layer deserves 
a mor e detailed discussion (see, e.g., iBeard fc WiescheTI 
(|2003l) ). Early calculations indicated that thermonuclear 
explosive burning produced mostly 56 Ni w hich then con- 
verted into 56 Fe by the electron captures (jTaam I . Il982l 
lAvasli fe Joss I . Il982t) . This picture had to be revised in the 
light of more recent simulations of nuclear evolution dur- 
ing cooling following the temperature peak of ~ 2 x 10 9 K. 
These simulations have shown that after a few minutes af- 
ter the initial temperature peak, nu clear ashes contain a 
mixture of nuclei with A ~ 50 - 100 (jSchatz et al. I . l200lh . 

Reactions in accreting crust. During accretion, the crust 
is a site of exothermic reactions in a plasma which is far 
from the catalyzed state. A layer of ashes of X-ray bursts is 
compressed under the weight of cumulating overlaying ac- 



creted layers. The ashes of density greater than 10 5 g cm -3 
are a plasma composed of nuclei immersed in an electron 
gas. The temperature in the deeper layer of thermonuclear 
ashes (a few 10 s K at the depth of a few meters) is too 
low for the thermonuclear fusion to proceed: the nuclei 
have Z ~ 30 — 50, and the fusion reactions are blocked 
by the Coulomb barriers. Therefore, the only nuclear pro- 
cesses are those induced by compression of matter. These 
processes are: electron captures, and neutron emissions and 
absorptions. Compression of ashes results in the increasing 
density of their layer (and its increasing depth below the 
stellar surface), and leads to the electron captures on nu- 
clei and the neutronization of the matter. To be specific, 
let u s consider a neutron star with M = 1.4 M Q (Fig. 
39 in IChamel fc Haensel I (|2008t) ). For densities exceeding 
~ 5x 10 11 g cm -3 (neutron drip point in accreted crust, at 
the depth ~ 300 m) electron captures are followed by the 
neutron emission from nuclei. Therefore, apart from being 
immersed in electron gas, nuclei become immersed in a neu- 
tron gas. At the density greater than 10 12 g cm -3 (depth 
greater than ~ 350 m) the value of Z becomes so low, 
and the energy of the quantum zero-point motion of nu- 
clei around the l attice sites so high, that the pycnonuclear 
fusion (see, e.g.. IShapiro fc Teukolskvl (|1983[ )) of neighbor- 
ing nuclei might be possible. This would lead to a further 
neutronization of the considered layer of accreted matter. 
With increasing depth and density, the element of matter 
under consideration becomes more and more neutronized, 
and the fraction of free neutrons outside nuclei in the total 
number of nucleons increases. At the density 10 14 g cm~ 3 
(depth ~ 1 km) the crust matter element dissolves into a 
homogeneous n—p—e plasma, containing only a few percent 
of protons. 

Nuclear composition of accreted crusts. The composi- 
tion of a fully accreted crust (all the crust, including its 
bottom layer, being obtained by compression of an initial 
shell with p = 10 8 g cm~ 3 ) is calculated , assuming a simple 
mode l of one-component plasma, as in (jHaensel fc Zdunild . 
I2003D . Two initial values of the mass number of nuclei in the 
X-ray bursts ashes are assumed, Ai = 56, 106. Thermal ef- 
fects are neglected. Extensive numerical simulations of the 
nuclear evolution of multi-component ashes, assuming large 
reaction network, and taki ng into accoun t temp erature ef- 
fects, were carried out by (jGuota et al. I . l2007t ). However, 
the latter calculations were restricted to densities less than 
10 11 g cm -3 , where the EOS is not significantly different 
from EOS [cat] , albeit the nuclear composition is very dif- 
ferent from the catalyzed matter one. 

A model of the EOS of a fully accreted crust will be 
denoted as EOS[acc..A], where A refers to the assumptions 
underlying the model: composition of X-ray bursts ashes, 
types of reactions included in the model, etc. 

2.3. Equations of state 

To disentangle the effect of formation scenario from the 
EOS of the crust, one has to compare the EOSs calculated 
not only for the same nuclear Hamiltonian, but also for the 
same model of nuclei and of nuclear matter in nuclei and 
neutro n matter matter ou t side nu clei. In this respect anal- 
ysis of (jHaensel fc Zdunik1 . ll990bl) was not correct, because 
in that paper the EOSs of accreted and catalyzed crusts 
were based on different nuclear models. In the present pa- 
per we use consistently the compressible liquid drop model 
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of nuclei of (jMackie &; Bavm 1 . 119771) : this model was used in 
all previous calculations of accreted crusts. We calculated 
several EOSs of accreted crust corresponding to different 
Ai of X-ray bursts ashes and for different models of pyc- 
nonuclear reaction rates. These EOSs are plotted in Fig.[TJ 
where we show also EOS [cat] for the same model of the 
nucleon component (nuclei and neutron gas) of dense mat- 
ter. For p < 10 11 g cm" 3 all EOSs are very similar, and 
therefore we display EOS plots only for p > 10 11 g cm -3 . 

Stiffness of the EOS. Significant differences in this prop- 
erty of the EOS exist in the inner crust, from the neutron 
drip point on, and up to 10 13 g cm" 3 . In this density range 
EOS [ace] are significantly stiffer than EOS [cat]. The dif- 
ference starts already at the neutron drip point, which in 
the accreted crust is found at a higher density. The soft- 
ening that follows after the neutron drip is much stronger 
in EOS [cat] than in EOS[acc]. Then, for p > 10 13 g cm -3 
EOS[acc] converge (from above) to EOS [cat], because at 
such high density the nuclei play a lesser role in the EOS 
and the crust pressure is mostly determined by the neutron 
gas outside nuclei. 

Jumps and smoothness. EOS[acc] shows numerous pro- 
nounced density jumps at constant pressure, to be con- 
trasted with smooth curve of EOSfcat]. Actually, both 
features are artifacts of the models used. For example, 
EOS [cat] is completely smooth, because we used a com- 
pressible liquid drop model without shell correction for neu- 
trons and protons, and treated neutron and proton num- 
bers within the Wigncr-Scitz cells as continuous variables, 
in which we minimized the enthalpy per nucleon at a given 
pressure. For EOS[acc] we used one-component plasma ap- 
proximation, with integer neutron and proton numbers in 
the Wigner-Seitz cells. As we neglected thermal effects, we 
got density jump at each threshold for electron capture. 
Had we used a multicomponcnt plasma model and included 
thermal effects and large reaction network, the jumps in the 
EOS[acc] would be smoothed fsee lGupta et al. I (|2007l) ). 

Pycnonucle ar fusion. The process of py cnonuclear fu- 
sion (see, e.g., (jShapiro fc TeukolskvL 11983^ may proceed 
after electron captures followed by neutron emission - a re- 
action chain which results in decreasing Z and A. In our 
simulations, based on the one-component plasma model, 
pycnonuclear fusion proceeds as soon as the time r pyc 
(inverse of the pycnonuclear fusion rate) is smaller than 
the timescale r acc needed for a matter element to pass, 
due to accretion, across the (quasistationary) crust shell 
with specific (A,Z). This time can be estimated as r acc = 
Mshoi\{A, Z)/M. Typical values of the mass of shells in the 
inner accreted crust are M s hoii ~ 10" 5 M Q . At accretion 
rates M ~ 10~ 10 — 1 0~ 9 M R yr -1 formulae used to calc ulate 
the fusion rates in (jHaensel fc Zdunikl . Il990al l2003ft lead 
to pycnonuclear fusions proceeding at p ^ 10 12 g cm" 3 . 

is plagued by huge 



However, theoretic al evaluation of r pvc 
uncertainties (see lYakovlev et al. I ( 2006[ )). It is not cer- 
tain whether pycnonuclear fusions do indeed occur below 
10 13 g cm -3 . If they do not, then EOS[acc] is quite well 
represented by models with pycnonuclear fusion switched 
off. The two extremes - pycnonuclear fusions starting at 
10 12 g cm -3 , and pycnonuclear fusion shifted to densities 
above 10 13 g cm -3 , correspond to curves EOS[acc.Ai.P] 
and EOS[acc.Ai.NP] in Fig.[TJ 

Bottom crust layer with p > 10 14 g cm" 3 . Our plots of 
EOS in Fig. [TJ are restricted to p < 10 14 g cm -3 . For the 



densities above 10 14 g cm" 3 and up to the crust bottom 
density pt>, the precision of our models lowers significantly 
compared to the precision of our EOS at lower densities. 
Fortunately, the contribution of the bottom layer of the 
crust (10 14 g cm" 3 < p < Pb) to the difference i? acc — i? cat 
is negligible and therefore the uncertainties in the crust 
EOS at highest densities do not affect our main results. 
This favorable situation will be justified in Sect. 3.2. 
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Fig. 1. (color online) EOS for the crusts calculated us- 
ing the Mackie-Baym model of nuclei and neutron mat- 
ter. Explanation of the labels: cat - cold catalyzed matter; 
acc. 56. P - initial X-ray burst ashes A\ = 56, pycnonuclear 
fusion switched on; acc.56.NP - initial X-ray burst ashes 

M 

For further details see Sect.l 



56, pycnonuclear fusion switched off till 10 13 g cm 3 . 



3. Effect of the change of the crust EOS on the 
star radius 

3.1. Accreted vs. catalyzed crust: numerical results for 
non-rotating models 

We matched EOSs of the crust described in the preced- 
ing section to several EOS of the liquid core. We checked 
that the difference AR(M) = R acc (M) ~ R cat (M) docs not 
depend on the details of matching of the crust-core inter- 
face. On the other hand, AR(M) can be shown to depend 
on the EOS of the core via the factor (R c &t/r g — l)i? ca t, 
where r g = 2GM/c 2 = 2.96 M/M Q km is the gravitational 
(Schwarzschild) radius. These two properties will be derived 
using the equations of hydrostatic equilibrium in Sect. 13.21 
In Fig. [2] we show the M(R) relation for neutron stars 
with catalyzed and accreted crust (EO S[acc.56.P]). For the 
liquid core we use the SLy EOS of (iDouchin fc HaenselL 
l200lt). For M = 1.4 M Q we get AR « 100 m. The value of 
Ai? grows to 200 m if M decreases to 1 M Q . It decreases 
to 40 m for 1.8 M . 
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Fig. 2. (Color online) Solid lines give gravitational mass 
M versus radius R for non-rotating neutron stars with 
accreted [acc.56.P] and catalyzed crust. Dashed red lines 
correspond to the approximation of the crust thickness by 
Eq. ([5]). Zoomed region illustrates a very high accuracy of 
this approximation. For the explanation of the procedure 
of getting dashed lines see the last fragment of Sect. 3.2. 



3.2. An analytic approximation for AR 

Let us consider hydrostatic equilibrium of a non-rotating 
neutron star. Let the (circumferential) radius of the star 
be R and its (gravitational) mass be M. The pressure at 
the bottom of the crust is P\, . Let the mass of the crust be 
M CT , We assume that M CI /M <C 1, but we account for the 
radial extension of the crust and the radial dependence of 
the pressure and of the d ensity within the cru st, following 
the method formulated in (jZdunik et al. I . l2008f) . We use the 
fact that within the crust P/c 2 < p and Airr^P/Mc 2 < 1. 
We define a dimcnsionless function of pressure within the 
crust (0 < P < P h ) 



X(P) = 



dP' 



p{P')c 2 



(1) 



The function x(P) is determined solely by the EOS of the 
crust. 

Using the Tolman- Oppenheimer- Volkoff e q uation of hy- 
drostatic equilibrium (Sh apiro fc Teukolskvi (11983D ). and 
neglecting within the crust P/c 2 compared to p and to 
M/(47rr 3 ), we can go over in Eq. ([T]) to the radial coor- 
dinate r, getting 



X [P(r)} 



1, 



-In 



rJR 



1 - r g /r 



(2) 



where r g = 2GM/c 2 . 

The dimensionless function % <C 1 can be treated as 
a small parameter in systematic expansions of the crust 
thickness. This function increases monotonously with P, 



from zero to ~ 10~ 2 (upper panel of Fig. [3]). Therefore, an 
expression for the crust layer thickness, from the surface 
to the pressure at the layer bottom P, i(P) = R — r(P), 
obtained from Eq. @ in the linear approximation in x- is 
expected to be very precise, 



t(P) = 2 X (P) 



R 



1 R 



(3) 



Our aim is to evaluate the change in the radius of neu- 
tron star of mass M and radius R, when EOS [cat] is re- 
placed by EOS[acc]. These EOSs differ for pressures below 
Pi = 10 32 erg cm- 3 (Sect.^). To calculate P acc -P cat it 
is therefore sufficient to know the difference in values of x 
at P\. We introduce a function A%(P), measuring the dif- 
ference between two EOSs of the crust for pressures from 
zero to P, 



A X (P) =Xacc(P)-Xcat(P) 



(4) 



In Fig. [3] we show functions x{P) an d A%(P) for several 
EOS of accreted crust. In spite of the jumps in density at 
constant pressure, characteristic of EOS of accreted crusts 
(Fig.[l|, both x an d Ax are smooth functions of P. This is 
due to the integration over P' < P in the definition of x(P)- 
As seen in Fig. [3j dependence of x and Ax on the particu- 
lar scenario underlying EOS[acc] is very weak. Additionally, 
Ax is nearly constant for P > 2 x I0 31 erg cm -3 , be- 
cause EOS [cat] and EOS[acc] converge at high pressures 
(see Fig. [TJ . Therefore, Ax(Pi) is a very good approxima- 
tion for Ax(-Pb)- 

Our final formula, obtained using Eqs. ((3j) , , combined 
with approximations explained above, is 



R* 



R cat ~2A[ X }-(R/r s -l)R, 



(5) 



where A[x] = Ax (Pi) and R = P cat . 

The precision of the approximation ((5|) is illustrated in 
Fig. [2 The dashed line for the accreted crust was obtained 
from exact P ca t(M) using Eq. ([3]) with R = R ca t- On the 
other hand, the dashed line for catalyzed crust was obtained 
from exact R acc (M) using Eq. ([5]) with R — P acc . 

Up to this point, we neglected contribution of elas- 
tic strain to the stress tensor within the crust, and used 
the perfect liquid approximation in the equations of hy- 
drostatic equilibrium. Recent numerical simulations indi- 
cate, that on the timescales of years, and longer (which 
are of astrophysical interest), the breaking stress of the 
i on crystal of the crust is ~ I0~ 3 of the crust pressure 
(jChugunov fc Horowitz I . l2010f) . This is much smaller than 
the difference between EOS[acc] and EOS [cat]. Therefore, 
as long as we restrict ourselves to timcscalc longer than a 
few years, contribution of clastic strain to AR(M) can be 
neglected. 

4. Combined effects of accretion and rotation on 

AR(M) 

Millisecond pulsars are thought to be recycled old 
("dead") pulsars, spun- up by ac cretion of matter 
from their companion in LMX Bs (lAlpar et all . If98l 
iBhattacharva fc van den Heuvel I . Il991h . This scenario is 
corroborated by discovery of rapid X-ray pulsations with 
frequencies up to 619 Hz, in more than a dozen LMXBs. 
The most rapid millisecond pulsar (isolated one) rotates at 
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Fig. 3. (color online) Upper panel: functions x(P) f° r 
EOS[acc.56.P] (dashed, red), EOS[acc.l06.NP] (dotted, 
blue) and EOS [cat] (solid, black). Lower panel: func- 
tion Ax(P) (Eq. |U for EOS[acc.56.P] (dashed, red) and 
EOS[acc.l06.NP] (dotted, blue) 



716 Hz. To reach such high rotation frequencies, neutron 
star had to accrete some five times the mass of the crust of 
a 1.4 Mq neutron star. Therefore, one has to conclude that 
millisecond pulsars have fully accreted crusts. 

In the case of rotating stars we will define i? a cc — -Rcat as 
the difference in circumferential equatorial radii. Rotation 
increases the difference i? acc — -Rcat compared to the static 
case. A rough Newtonian argument relies on the propor- 
tionality of the centrifugal force to the distance from the 
rotation axis. Therefore, centrifugal force that acts against 
gravity is largest at the equator. In order to get precise value 
of the effect of rotation on i? acc — i? ca t, we performed 2-D 
simulations of stationary configurations of rigidly rotating 
neutron stars with catalyzed and accreted crusts. The cal- 
culations have been performed using the rotstar code from 
the LORENE library (http : //www. lorene . obspm.fr Our 
results, obtained for / = 716 Hz (maximum frequency mea- 
sured for a pulsar) and for a half of this frequency, are shown 
in Fig. gl 

Consider hydrostatic equilibrium of the crust of rotat- 
ing star in the equatorial plane. In the Newtonian approx- 
imation, the ratio of centrifugal force to the gravitational 
pull at the equator is, neglecting rotational deformation, 
Vl 2 R 3 /GM, where fl is angular frequency of rigid rotation. 
Notice that this expression is exact in the quadratic ap- 
proximation in because increase of R due to to rota- 
tion introduces higher powers of f2. We propose modeling 
centrifugal-force effect within crust by modifying Tolman- 
Oppcnheimer-Volkoff equation in the equatorial plane in 
following way: 



dP 
dr 



GMp 



1 - 



2GM\ 



(1 - an 2 R 3 /GM) (6) 



where M is the total mass of the star and a is a numerical 
coefficient to be determined by fitting the exact results of 
2-D calculations. The validity of Eq. ^ relies on the small- 
ness of the rotational flattening of the liquid core compared 
to that of crust. We included terms quadratic in fl and 
we used standard approximation valid for a low-mass crust 
and applied in the preceding section. Equation © could be 
solved explicitly assuming constant M and taking into ac- 
count the changes of r throughout t he crust, as in the case of 
non-rotating star in (IZdunild . I2002T ) . However, to be consis- 
tent with the approach leading to formula ([3]) , we prefer to 
use a solution method based on the smallness of \- Rotation 
effect enters via a constant factor (l — afl 2 R 3 /GM). This 
results in a simple relation between 51 = and ft > dif- 
ference in radii: 



Ai? oq (f2) — i?eq,acc(^) ~ #eq,cat(^) 

_ 2A\x\-(R/r s -l)R 



Ai?(Q = 0) 



l-aQ 2 R 3 /GM l-an 2 R 3 /GM 



(7) 



The dimensionless parameter a has been determined nu- 
merically by fitting formula ([7]) to the exact 2-D results ob- 
tained in General Relativity using LORENE numerical li- 
brary. We found that a ~ 4/3. As one sees in Fig. [H our ap- 
proximate formula for Ai? cq (f2) works extremely well. The 
actual f2-dependence of Ai? eq is stronger than quadratic 
because ft 2 appears in the denominator in Eq. (|7|). 
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Fig. 4. Dependence of Ai?(M) = i? acc (M) - i? cat (M) on 
the frequency of rotation of neutron star, /. Equations of 
state as in Fig. [5] Upper curve - / = 716 Hz. Middle curve 
- / = 358 Hz. Bottom curve - non-rotating. 



5. Summary and conclusion 

In the present note we studied the effect of the formation 
scenario on the mass-radius relation for neutron stars. For 
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a given M, a star with an accreted crust has larger radius, 
by AR(M) , than a star built of catalyzed matter formed in 
stellar core collapse. We derived an approximate but very 
precise formula for AR(M), valid for slowly rotating neu- 
tron stars. AR(M) factorizes into a prefactor depending 
solely on the EOS of neutron star crusts formed in differ- 
ent scenarios and a simple function of M for a given EOS 
of the core. We studied the dependence of the difference 
between the equatorial radii on the angular rotation fre- 
quency, Ai?eq(0). We derived an approximate formula for 
A-Req(Q), that reproduces with high precision Ai? eq even 
for neutron stars rotating at 716 Hz, highest rotation fre- 
quency measured for a radio pulsar. 

We found that an accreted crust makes the radius of 
a 2 M Q — 1 M Q star some 50 — 200 m larger than in the 
standard catalyzed matter case. Highest hopes of a simul- 
taneous measurement of a neutron star M and R are, in 
this deca de, associated with high resolution X-ray spec - 
troscopy (jArzoumian et all . 120091 iPaerels et al. I . l2009f) . 
Unfortunately, expected uncertainty in determining R(M) 
is ±5%. It significantly exceeds effects of formation scenar- 
ios, calculated in the present note. 
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